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Abstract 

Obtaining accurate Channel State Information (CSI) at the transmitters (TX) is critical to many 
cooperation schemes such as Network MIMO, Interference Alignment etc. Practical CSI feedback and 
limited backhaul-based sharing inevitably creates degradations of CSI which are specific to each TX, 
giving rise to a distributed form of CSI. In the Distributed CSI (D-CSI) broadcast channel setting, 
the various TXs design elements of the precoder based on their individual estimates of the global 
multiuser channel matrix, which intuitively degrades performance when compared with the commonly 
used centralized CSI assumption. This paper tackles this challenging scenario and presents a first analysis 
of the rate performance for the distributed CSI multi-TX broadcast channel setting, in the large number 
of antenna regime. Using Random Matrix Theory (RMT) tools, we derive deterministic equivalents 
of the Signal to Interference plus Noise Ratio (SINR) for the popular regularized Zero-Forcing (ZF) 
precoder, allowing to unveil the price of distributedness for such cooperation methods. 

Index Terms 

Multiuser channels. Cooperative communication, MIMO, Feedback Communications 

I. Introduction 

Network (or Multi-cell) MIMO methods, whereby multiple interfering TXs share user mes¬ 
sages and allow for joint precoding, are currently considered for next generation wireless net¬ 
works [1]. With perfect message and CSI sharing, the different TXs can be seen as a unique 
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virtual multiple-antenna array serving all RXs in a multiple-antenna broadcast channel (BC) 
fashion, and well known precoding algorithms from the literature can be used [2]. Joint precoding 
however requires the feedback of an accurate multi-user CSI to each TX in order to achieve 
near optimal sum rate performance [3]. 

Although the case of imperfect, noisy or delayed, CSI has been considered in past work [3], 

[4] , literature typically assumes centralized CSIT, i.e., that the precoding is done on the basis of 
a single imperfect channel estimate which is common at every TX. Although meaningful in the 
case of a broadcast with a single transmit device, this assumption can be challenged when the 
joint precoding is carried out across distant TXs linked by heterogeneous and imperfect backhaul 
links or having to communicate without backhaul (over the air) among each other, as in the case 
of direct device-to-device cooperation. In these cases, it is expected that the CSI exchange will 
introduce further delay and quantization noise. It is thus practically relevant for joint precoding 
across distant TXs to consider a CSI setting where each TX receives its own channel estimate. 
This setting is referred to as distributed CSI (D-CSI) in the rest of this paper. 

From an information theoretic perspective, the study of transmitter cooperation in the D-CSI 
broadcast channel setting raises several intriguing and challenging questions. 

First, the capacity region of the broadcast channel under a general D-CSI setting is unknown. In 

[5] , a rate characterization at high SNR is carried out using DoF analysis for the two transmitters 
scenario. This study highlighted the severe penalty caused by the lack of a consistent CSI shared 
by the cooperating TXs from a DoF point of view, when using a conventional precoder. It was 
also shown that classical (regularized) robust precoders [6] do not restore the DoF. Although a 
new DoF-restoring decentralized precoding strategy was presented in [5] for the two TXs case, 
the general case of more than 2 TXs remains open. At finite SNR, the problem of designing 
precoders that optimally tackle the D-CSI setting is open for any number of TXs. The use of 
conventional linear precoders that are unaware of the D-CSI structure is expected to yield a loss 
with respect to a centralized (even imperfect) CSI setting. Yet the quantifying of this loss in the 
finite SNR region has not been addressed previously. This is precisely the question addressed 
by this paper. 

We study the average rate achieved when the number of transmit antennas and the number of 
receive antennas jointly grow large with a fixed ratio, thus allowing to use efficient tools from 
the field of RMT. Although RMT has been applied in many works to the analysis of wireless 
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communications [See [7]-[ll] among others], its role in helping to analyze cooperative systems 
with distributed information has not been highlighted before. 

Our main contribution consists in providing a deterministic equivalent for the average rate 
per user in a D-CSI setting where each TX receives its own estimate of the global multi-user 
channel matrix with the quality (in a statistical sense) of this estimate varying from TX to TX. 


II. System Model 


A. Transmission Model 


We study a so-called Network MIMO transmission where n TXs serve jointly K receivers 
(RXs). In order to simplify our analysis we restrict ourselves to linear precoding structures. Each 
TX is equipped with Mtx antennas such that the total number of transmit antennas is denoted by 
M = nM^x- Every RX is equipped with a single-antenna. We assume that the ratio of transmit 
antennas with respect to the number of users is fixed and given by 




( 1 ) 


We further assume that the RXs have perfect CSI so as to focus on the imperfect CSI feedback 
and exchange among the TXs. We consider that the RXs treat interference as noise. The channel 
from the n TXs to the K RXs is represented by the multi-user channel matrix H G . 
The transmission is then described as 


( 2 ) 


where y* G C is the signal received at the i-th RX, hf = efH G is the channel from 

all transmit antennas to RX i, and rj = [r]i,... G is the normalized Gaussian noise 
with its elements i.i.d. as CJ\f{0, 1). 

The transmitted multi-user signal x G is obtained from the symbol vector s = ..., 


yi 


’/l^X 




= Hx + T] = 


+ 


Vk 




Vk 


K\ ^ 


.Kxl 


with its elements i.i.d. Mc{0, 1) as 


X = Ts = 


LR 


Si 


Sr 


( 3 ) 
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with T G being the multi-user precoder and = Te* G being the beamforming 

vector used to transmit to RX i. We consider for clarity the sum power constraint ||T||| = P. 
Our main figure-of-merit is the average rate per user 

1 ^ 

R^-5^E[log2(l + SINR,)] (4) 

k=l 


where SINR^ denotes the Signal-to-Interference and Noise Ratio (SINR) at RX k and is defined 
as 


SINRa. = 


\hft 


k ^k\ 


1 + Y.f=l/^k 


(5) 


B. Distributed CSIT Model 

In the distributed CSIT model studied here, each TX receives its own CSI based on which 
it designs its transmission parameters without any additional communication to the other TXs. 
Specifically, TX j receives the multi-user channel estimate G and designs its transmit 

coefficients G c^^xxi solely as a function of For ease of exposition, we assume that 
the imperfect multi-user channel estimate is modeled by 

- (a(f))2H + (6) 

with G having its elements i.i.d. Afc{0, 1). We then denote by the kth row 

of i.e., the estimate at TX j of the channel from all the transmit antennas to RX k. The 
approach described in this work extends to a more general information structure allowing for a 
non uniform description quality at each TX as well as correlation in the multi-user channel. The 
calculations will be provided in an upcoming work. 

Remark 1. The D-CSI model encompasses the imperfect centralized CSI model by taking n = 

1 . □ 


C. Regularized Zero Forcing with Distributed CSIT 


We address the performance of a classical MISO broadcast precoder, namely regularized ZF 
[6], [12], when faced with distributed CSIT in the large system regime. Hence, the precoder 
designed at TX j is assumed to take the form 





A + MalM 


(7) 
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with regularization factor a > 0. We also define 

M 

such that the precoder at TX j can be rewritten as 


■v/'l/O) 


The scalar corresponds to the power normalization at TX j. Hence, it holds that 

^ 0 ) = + MalMY^ 


Upon concatenation of all TX’s precoding vectors, the effective global precoder denoted by 
is equal to 

■*^1 -*-rZF 

TpH'Ti{2) 
rpDCSI A 2 -^rZF 
-■-rZF ~ 

■rH rr\{K) 

where E? G c^^xxm defined as 




*AfxX ^ 0 — 1)-^TX ^M'yx. ^ — j)^TX. 


We furthermore denote the kth column of (used to serve RX k) by 

Although the finite SNR rate analysis under the precoding structure (11) and the distributed 
CSI model in (6) is challenging in the general case because of the dependency of one user 
performance on a all channel estimates, some useful insights can be obtained in the large antenna 
regime, as shown below. 


III. Deterministic Equivalent of the SINR 

The Stieltjes transform has proven very useful many times in the analysis of wireless networks 
[9], [10] and we will also follow this approach. Hence, our approach will be based on the 
following fundamental result. 

Theorem 1. flO], [131 Consider the resolvent matrix Q = f with the matrix H 


Theorem 1. [10], [13] Consider the resolvent matrix Q = 1^- 
defined according to Section 11 and a > 0. Then the equation 


T=^tr + 
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admits a unique fixed point which we will denote by 6 in the following. Let 

and let the matrix U be any matrix with bounded spectral norm. Then, 

itr(UQ)-itr(UQ„)--^0. (15) 

Aj Aj K,M^oo 

The fixed point 5 can easily be obtained by an iterative algorithm given in [9], [11] and 
recalled in Appendix A for the sake of completeness. Using this theorem and the definition of 5, 
we can now state our main result. 


Theorem 2. Considering the D-CSI BC described in Section II, then SINR^. — SINR^ —)■ 0 with 
SINR^ defined as 


SINR2 ^ 


with /fc G R given by 


4° = Z 


^ yi - 
4” + ¥ 


n + 26 (—1 + n + 


^ (1 + 6W 

j=i 

+5' (-1 + n + 

+ (a(^'))yi - (a(^'))y(a0’))2 - 


r° 6 


n n 

(1 + 6Yn 

+(i (—1 + 


(i)U ^ 


-2 + {a'^^>Y+ {a 


while r° G R and V° -i G R are respectively defined as 


52 


■^O A /3(l+'5) 




■po ^ 

^ jj' ~ 


(1 - ^) + 


1 _ i. 


/? (l+5f. 

Y (l-((T0'))2)(l-(cr0''))2)52 

^(T+5) 


(1-5) + 


52 

(1+5) 


1 _ 1 (1 _ (^(.)) 2 ) (1 _ (^(/)) 2 ) 


52 


(1+5)^ 


(16) 


(17) 


(18) 


(19) 
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IV. Proof of Theorem 2 

Our calculation is built upon results from both [9] and [10]. We also make extensive use of 
elassieal RMT lemmas reealled in Appendix A. Note that Lemma 7 and Lemma 8 are novel. In 
partieular, Lemma 7 extends [10, Lemma 15] and is an interesting result in itself. 

During the ealeulation we use the notation a; x u to denote that x — y ———)■ 0. 

K,M^oo 

A. Deterministic Equivalent for 

We start by finding a deterministie equivalent for In faet, a deterministic equivalent for 
is provided in [9]. However, it ean also be obtained using Lemma 7 with = 0, 

whieh gives 

X T°. (20) 

Looking at the definition of r° in (18), it ean be noted that, as expeeted, this deterministie 
equivalent does not depend on 


,H+DCSI . 


B. Deterministic Equivalent for k- 

Turning to the desired signal at RX k, we ean write 


lH.DCSI 
'‘'fc LZF.fc 


where we have defined 


1 ^/P 
^ M 






(a) 


(b) 


(c) 


id) 



j V" J ■■'k 

1 £, (i) 




7=1 


-E 


[O' 


ij))‘ 


[k]> "'k 

-n. 


M 





7 = 1 


E 

7 = 1 


1 


[k]. 


((T 


(7)): 


Mtr (E,EH(C 


(7)'|-l 

[k]) 




5 


7 = 1 


1 + 5 


r^ij) A 

^[k] - 


Vj 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


M 


+ ciIm, 


(26) 



with 


fH 


(i)\H A 
[k]) - 


h\ 


U) 


lAj) iAj) 
"'k-l "'fc+1 


h 


Vj. 


(27) 


Equality (a) follows then from Lemma 2 and the use of the deterministie equivalent derived 
for ip) from Lemma 4, (c) from Lemma 3, (d) from Lemma 5 and the fundamental 
Theorem 1. Note that 5 is defined in Theorem 1. It follows then directly that 


i=i 


(a0))2 


^2 


(l + <5) 


2 • 


(28) 


C. Deterministic Equivalent for the Interference Term 

Our first step is to write explicitly the interference term using the definition of in (18) 

and replace by its deterministic equivalent. 


K 


''rZF,^ I 


e=l,£jtk 

.Ht-iDCSI/t-iDCSInH 


_ ^H.aDCSW.aDCSI\H^ 

— '''fc -'-rZF l-'-rZF ) '^k '‘'fc‘'rZF,fcl‘'rZF,fcJ '^k 


^ n n p 


M2 ^ 


j= 

p 

1 

TO 

M2 

P 

1 

r° 

M2 




ywVW) 


yy h^E,E«(c“)-'(Hg)«H«(Cb'')-'%E'!h> 

jr = l j'=l 


(29) 

(30) 

( 31 ) 

(32) 


((C(^))-i - (Cg)-i) (H[g)«H[^;)(C(^'))-iE,.E»/i,. (33) 


'[k]. 

j=i j'=i 

To obtain a deterministic equivalent for the second summation in (33) we use the following 
relation 


(c«')-‘ - (cg)-‘ = (€«>)-■ (c|)’ - ca>) (c“)-‘ 




(34) 


( 35 ) 
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Inserting (35) in (33) yields 
ik 

n 1 n n 


p 

1 

f° 

M2 

p 

1 

r° 

M3 

p 

1 

r° 

M3 

p 

1 

r° 

M3 

p 

1 

r° 

M3 


j = l j'=l 
n n 


j=i j'=i 
n n 




y; y; h«E,E«(C“)-' (C®)-'(H|'>)«H«(Cb'')-‘E,.E«ft, 

j = l j' = l 


y^y^h?EyE«(c«)-' 

j=i ji=i 


Sfcfht 


(c;^;)-'(H|',')«Ho;>(ce'>)-'E,,E«hi 




Z Z h?EiE»(Cl'>)-' [fe*cy(4'>)«] (C“)-'(H“)«Hg>(C<-’'>)-'E,,E»h, 

j = l j' = l 


= A- B -C - D - E. 


(36) 

(37) 


We proceed by calculating each of the 5 terms in (37) successively, using in particular Lemma 8: 


A = 


P 1 

f°ii^ 


Y Y '»?E,EH(cW)-'(H“)«H«;'(C'^'l)-‘E,.E«h, 


(38) 


j = l j'=l 
n 


tr (e,,e«e,e«(c®)-hh“)«P|„h“|'(c;,| 




- c, 


j=i j'=i 

(.,)tr (EjL, vv.[fc] 


M2 


,^tr (E,EH(Cg)-HHg)HP[,,Hg)(C(^;VM tr (E,,EH(cgVM 1 + 


,.tr((c[^3) 1 


+ (4^ 0 


tr E,EH(C ,, 

(y) 2 V ^ w- 


M2 

(HS)"PwH5’(Cj,| 


M 


1 + 




M 


tr E,vEH(C 






/ M 


M2 


M 


1 + 


-MM' 


M 


(39) 


From the unitary invariance of the distribution of H and it can be shown that 
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where the last equality follows directly from applying Lemma 7. Inserting (41) in (39) and using 
the fundamental Theorem 1 yields 


-yy- 

r° ^ ^ n 

j=i j'=i 


1 . . _ ^(T) hi' 0 -I- ~r <^1 0 j,j' ^ ^ 

° n n 1 + 6 ^ n nl + 6 


We then proceed similarly for the remaining 4 terms: 


r° M3 


n n 


j = l j'=l 
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c = 


TD 1 

j = l j'=l 


Si) 


P 1 




i=i 


Si) 


(4'’)"(C“)-‘(H[i,’)’‘H|«(C<-’'>)-'E,,E>, 

(46) 


n 


i=i 


(,.)“'(EiE«(cg)-‘) Kisfr: 


M 


1 + 


tr((C«|)-l 

M 


Si) Si) 

Cl C2 


,jtr ((CW)->(h“)«P|„h“(CW)->) tr (e,E«(C«)-‘) 


M2 


M 


1 + 




M 


— C. 


(i) 


tr ((CW)->(Hg)«P|r,H“(cW)->) tr (e,E 4C«)-4 , ^ „a)llMM 


l + c“ 


M 


M2 


M 


1 + 


tr((C«})-l 
M 


P ^{j)S 6 f Si) Si)-po^ ^ (i)po1 + *21 

nl + (5 r' ^ nl + (5 n 1 + 5 


5^4^ (-1)( 


1=1 


(47) 


(48) 


(49) 


D = 


P 1 
f°M3 


n n 

EE [+E,E«(C<-’I)-‘«“] . 4'. [hJ(C|'+'(H;4++''>(Cl«)-'E,.E»h; 


1=1 l'=l 


jj n n 

^EE' 


.( 1 ) 


-1)5 


( 1 ) 

2 


tr(E,EH(c“)->) !+|! 


M 


1=1 l'=i 

rtr (E,,E“(C“')-‘(H«;)“PHHg'(C;,, 


1 + 




M 


M2 


,^tr ((cgi+HHi^i+Pw+Jlcg')-') tr (e,.EH(c|«)->) i + 


Si') 


41') N 2 


M2 


M 


1 + 


to(C-i) 


M 


tr ((C®)->(Hli+PwH|j;’(C«;>)->) tr (e,,E>1(C 


H/r4(l')^-l 


/ M 


n n 

^EE4’( 

1=1 i'=i 


M2 


-DcO''* 


M 


1 + 


tr 




M 


po 

1 , 1 ' 


n 1 + 5 \ n 


,0')ro 'll+ 4 '',! . fJf))2po s s 

“ 1 + ^ +(<^2 (Ir/-- 


n 1+5 


(50) 


• (51) 
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E = 


P 1 

f°ii^ 




(4'’)"(Cw')-‘(H!S)*‘H«(C«'>)-‘E,-E“h, 


j = l j' = l 


P 1 
f°M3- 




(52) 


(53) 


i=i 


E 

j=i 


„(i) 

*-2 


tr I E,E«(C«)-M 1 + 4- 


Si) 


M 


M 


1 + 




M 


Cl ^2 


tr ((C^l-HHgl-PwHglC^)-') tr (e,E«(cW)->) HlSfH: 


.(i) 


M2 


tr((c[^;)-i(H[^;)HP[,]H[^;(c 


[fc]. 




M 


riH/r^(i)\-i 


tr(E,E^.(C[,]j l + c)- ^ 


1 + 
P) 


tr 




M 


tr 


((gw)-’] 


M2 


M 


1 + 


tr((C«)-l 

M 




^ ^ n 1 + 5 
i=i 


(54) 


(55) 


nl + 5 1 + 5 

The final expression is obtained after inserting all the deterministie equivalents derived inside the 
interferenee expression (37). The eompaet expression of the theorem is obtained after algebraie 
manipulations using the software Mathematiea. 

Remark 2. It is important to differentiate the eases j = j' and j ^ j' when eomputing r°^,. 
Indeed, in the ease j = j', it holds 

TP, = r°. (56) 


□ 


V. Simulation Results 

We now verify using Monte-Carlo simulations the aeeuraey of the asymptotie expression 
derived in Theorem 2. We eonsider a network eonsisting of n = 3 TXs with a sum power 
eonstraint given by P = 10 dB and a = 1/P. We foeus in this work on the ease of ((j*^ 20)2 _ 
0.1, Vj = 1,..., n so as to emphasize the price of distributedness. 

In Fig. 1, we show the rate per user as a funetion of the number of users for a square 
setting where M = nMxx = K (i.e., /9 = 1) in the distributed CSIT eonfiguration where 
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Fig. 1: Average rate per user as a funetion of the number of users K with = 0.1, Vj. 


((j0'))2 = 0.1, Vj = 1,... , n. For eomparison purpose, we also show the rate per user obtained 
in the ease of imperfeet eentralized CSIT with = 0.1 and with perfect CSIT (i.e. 

(crCcsi )2 = Q or equivalently = 0,Vj). As note earlier, a deterministic equivalent for the 

centralized case is obtained using n = 1 in Theorem 2. 

The large system deterministic equivalents are shown to be useful with just 20 to 30 users 
and antennas. In addition, the cost of having distributed information is also highlighted by the 
losses compared to the centralized configuration. This shows the necessity to take properly 
into account the CSI configuration when designing the feedback scheme and the precoder. 
Considering the same network configuration, we show in Fig. 2 the average rate per user in 
terms of the regularization factor a. Quite interestingly it appears that the optimal regularization 
















14 



Fig. 2: Average rate per user as a function of a for = 0.1, Vj. 


factors (represented with the * marker) are not the same in the centralized and distributed CSI 
settings. 


VI. Conclusion 

We have studied in this work the joint transmission using regularized ZF in a distributed CSI 
configuration. Using RMT tools, an analytical expression has been derived to approximate the 
average rate per user. This expression becomes asymptotically exact in the large system limit 
where the number of transmit antennas and the number of receive antennas go to infinity at the 
same pace. This new deterministic equivalent reveals the cost related to not just CSI feedback 
limitation, but also backhaul sharing limitations and can be helpful to design more robust systems. 
The extensions to more general channel and CSI models are challenging and subject to ongoing 
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work. Note that the price of distributedness is evaluated here for a conventional precoder. This 
further motivates the development of novel precoding schemes that are more suitable to the 
distributed CSI setting. 


Appendix 


A. Classical Lemmas from the Literature 

Lemma 1 (Adapted from [9], [11]). Let a > 0 and {Sk)k>o be the sequence defined as 


5n = i 


Then Sh 


4 = ig tr (alM + for k > I 


S, with S being by construction a fixed point of (57). 


Lemma 2 (Resolvent Identities [10], [11]). Given any matrix H e qKxm^ 
kth row and G denote the matrix obtained after removing the kth row from 

H. The resolvent matrices of H and are denoted by Q = (H^H + q;Ia^) ^ and Qk — 
+ q;Ia^) \ with a > 0, respectively. It then holds that 

^ Ml + ^h^Qkhk 




1 + Ifb^Qkhk 


Lemma 3 ( [10], [11]). Let (A 7 v) 7 v>i, Aat G be a sequence of matrices such that 

limsup IIAatII < oo, and (x7v)Ar>i, xat G be a sequence of random vectors of Lid. entries 
of zero mean, unit variance, and finite eighth order moment independent of A tv- Then, 

^xj^AjvXAT - ^ tr (Aat) 0. (60) 

iV iV N^oo 

Lemma 4 ( [10], [11]). Let A^- G be a sequence of matrices such that 

limsup IIAatII < oo, and be random, mutually independent with i.i.d. entries of zero 

mean, unit variance, finite eighth order moment, and independent of A at. Then, 

^x^AatYat > 0. (61) 

iV N^OO 
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Lemma 5 ( [9], [11]). Let Q and be as given in Lemma 2. Then, for any matrix A, we have 

tr (A (Q - Qfc)) < ||A|| 2 . (62) 

Lemma 6 ( [9], [11]). Let U,V,© be of uniformly bounded spectral norm with respect to 
N and let V be invertible. Further, define x = &z and y = Qq where z,q G have 
i.i.d. complex entries of zero mean, variance 1/N and finite 8th order moment and be mutually 
independent as well as independent o/U, V. Define cq, ci, C 2 G R’*' such that cqCi — > 0, and 

let u ^ ^ tr (©V“^) and u' = ^tr (©UV“^). Then we have: 

HN-1„ u'{l + ciu) 


x^U (V + coxx^ + ciyy^ + C 2 xy^ + C 2 yx^) x - 


(cqCi - + (co + Ci)u + 1 


^ 0 


(63) 


as well as 


x^U (V + coxx^ + ciyy^ + caxy^ + cayx^) ^ y - 


—C2UU 


(cqCi - + (co + Ci)u + 1 


^ 0 


(64) 


B. New Lemmas 


Lemma 7. Let H' and H" be two imperfect multi-user channel estimates as described in 

-1 . \-1 


Section II. Let Q' = 


M 


Oil 


M 


and Q" = 


M 


-|- Oil 


M 


with a > 0. Let 


A G C 


MxM 


be of uniformly bounded spectral norm with respect to M. Then, 


1 


M2 


/9(l + 5) 


(1-5) + 


52 + 

(1 + ^) 


with Cq = 1 — (t'2, Cq = 1 — (j"2, Qq defined as in Theorem 1, and Yq defined as 


.12 J! A 




_A 5(1+^) 
in — — 


(1 - <^) + (i^ 


1 _ r" 

^ 13 ^ 0^0 

Note that in the case where A = 1m, the result simplifies to 


(65) 


( 66 ) 


(67) 


Proof: We start by defining 


Q 


/ A 


M 


+ Cal 


M 


-1 


( 68 ) 
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with 

Hf ^ ... h',_^ ... h’^] . (69) 

We then define similarly Q" and Let us start by writing the simple equality 



We will now calculate separately each of the term Zj. Starting with Zi gives 



(75) 

1 ^ 1 

e=i 

(76) 

(a) 1 A 1 /if Q"AQo^' 

- M yMi + jihpq,"h'i 

(77) 

{b) 1 a/cqCq A tr (Q"AQo) 

(78) 

~ 1 + If w") 

(79) 

(rf) 1 ycoc[)V^tr(A) 

(80) 

j3 1 + 5 


where equality (a) follows from Lemma 2, equality (6) from Lemma 3, equality (c) from 
Lemma 5, and equality (d) from the fundamental Theorem 1. The following calculations are 
very similar and the same lemmas are used in the same way such that we will omit to mention 
explicitly the lemmas used. 
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Turning to Z 3 gives 
1 


^8 = 


M3 


M3 


tr (AQoH'“H'Q'H'«H"Q") 


K 


^tr {hfQ'H'^H”Q"AQoK) 


e=i 

K 


1 ^ tr Q'H'Hh"Q"AQ„/i9 


( 81 ) 

(82) 

(83) 


(a) 


e=i 

K 


M3 


1 ^ tr (^f Q'H'Hh"Q''AQo/i') 1 ^ tr Q''AQ„/r9 

2 - 1 + ihf Q',« 


£=i 




= ^4 + ^5 


(1 + iif^TQ'eK) (1 + 

(84) 

(85) 


with equality (a) obtained using Lemma 2. We also split the calculation in two and start by 
calculating Z 4 as follows. 

_ ^ tr (/if Q'/i'/if Q''AQo/i9 

K 


K 


Za = 


1 ^ tr (h«QlH«H;'Q';AQ„h;) 


AP 

1 

IP 


1 + ish’^Qf,K 


E 

£=1 

tr (Q^Hf H"Q"AQo) ± j—j 

^ l + ^tr(Qf mZ^VCo 


1 + i,/if 


r=i 


igli’(Q"AQo) 


,// M 


O'-O 


A^tr(Q'H'HH"Q"AQ„) K 


M l + Atr(Q') 

^ <^V^^^]^tr(A) 


M^ ° 


^=1 l + igti-(Qf 

^ tr (QO tr (Q"AQo) 


^ 0^0 ~ 


1 + l^tr(Q') 


/S 


1 + 5 


( 86 ) 

(87) 

( 88 ) 
(89) 
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Finally, it remains to calculate as 

' " ^ (1 + i)" 

1 tr (fe?Q;H«H?Q;fe?fe;HQ;AQ>',) 

AA‘ ^ (1 + if 


K 1 


M ^ 0 ' 


1 ^ tr (/if Q;AQ„/i') 

{i + W 

tr(Q'HfHM)igtr(Q''AQf 


( 1 +^r 


— (Q^) j4 (Q^O 14 (Q^ aQo) 
(1 + sf 


1 ^ , ,/aP tr (Q^HfH"Qf ^tr(QoAQo) A tr (Q^AQo) 


(l + <5r 




^ 1 , „ i.tr(Q-HfH--Qf i5^^tr(A) 1 tr (A) 

(l + <5f /S (l + <5f ' 

Adding all the Zi gives 


tr (AQ'H'“H"Q") 


1 

(3 1 + 5 


v/^5^^tr(A) 
/3(l + 5) 


52 c'c''^tr(Q'HfH"Qf 52 


c'qc'^S^ \ tr(A) 


/9 1 + 5 /9(l + 5)" 

52(l + v/^^tr(Q'HfH"Q; 


f3{l + 5y M 


It remains then to calculate the case A = to conclude the calculation. In that case, we have 


r' +'52 

CqCqO 


Hence, 


1 - ° ° , ^tr (Q'H'«H"Q") 


4tr(Q'H'HH"Q")x^ 


/3(l + 5) 


(1-5) + 


(l + <5) 


(1 - ^) + (i^ 


1 _ _£o£o^ 

m+5f 


Inserting (99) inside (98) concludes the proof. 
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Lemma 8 . Let L,R, A G be of uniformly bounded spectral norm with respect to M 

and let A be invertible. Let x, y have i.i.d. complex entries of zero mean, finite variance and 
finite 8th order moment and be mutually independent as well as independent o/ L, R, A. Then 
we have: 


with 


x"LA-iRx 

M 

x«LA-^Ry 

M 


UhK — ColfLltR- 


1 + CiU 

1 + U 


MLR + C1C2MLMR 


U 


1 + U 


+ C 2 MLMR 


— C2 MlMR- 


U 


1 + U 

1 + CiU 

1 + U 


A = A + cqxx^ + Ciyy" + C 2 xy" + C 2 yx 


H 


H 


H 


with Co + Cl = 1 and cqCi — C 2 = 0 , and 


A 

u — 


tr(A 


^ tr(LA 
Ml = 


-I'l 


^ tr(A iR) 
Mr = 


M ’ M ’ M 

Proof: Focusing first on the first equality gives 

—x^LA ^Rx- x^LA ^Rx 

M M 

= —x^LA"^ (A - A) A-^Rx 

= -^x^LA-i (coxx« + ciyy“ + C2yx« + csxy^) A’^Rx 
(a) 1 , „ , „ -I . tr (A“^R) 

— (cqX^LA ^X + C 2 X^LA y) ^ 


^ tr(LA-iR) 
MlR — - 


M 


M 

(b) tr (LA“^) tr (A^^R) 1 + Cl 
X -Co 


M 


tr( A- 


M 


M 


M 


1 + 


tr( A-l 


+ C 2 


tr (la tr (a ^R) 


tr( A- 


M 


M 


M 


1 + 


tr( A-l 1 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 


M ' M 

where equality (a) is obtained from using Lemma 4 and Lemma 3 and equality ( 6 ) follows from 



21 


Lemma 6. Similarly, we turn to the second equality to write 


—x^LA-iRy - —x^LA-^Ry 
M ^ M ^ 

=—x^LA^i (A - A) A-^Ry 
M y ) y 

= - ^x^LA"^ (cqxx^ + ciyy^ + Csyx^ + Caxy^) A"^Ry 

(“) 1 / H 1 H 1 \ tr fA“^R) 

x-^(cix«LA-V + C2X^LA-1x) ^ > 


M 


(6) tr (la tr (A ^R) 

XC 1 C 2 - 


tr( A- 


M 


M 


M 


1 + 


tr A- 


C 2 


tr (la tr (A ^R) 1 + Ci 


tr( A- 


M 


M 


M 


1 + 


tr A- 


(105) 

(106) 

(107) 

(108) 
(109) 


M ^ ' M 

where equality (a) is obtained from using Lemma 4 and Lemma 3 and equality (6) follows from 
Lemma 6. ■ 
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